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(57) Abstract: The invention relates to a method for computationally determining optical properties of a channel waveguide, accord- 
ing to which, for an incident ray entering an entrance surface, the distribution of intensity over an emergence surface is determined 
by adjoined partial rays that possibly split up. From an incident partial ray striking the outer surface, a predetermined algorithm 
determines a reflected main ray and, in so far as it is necessary according to the ensuing course of events, determines a number of 
scattered rays of a higher order, which depict partial rays and are recursively traced further. Each partial ray is either a primary ray or 
a secondary ray. When the incident partial ray is a primary ray, the reflected main ray is a primary ray once again, and the scattered 
rays are secondary rays. During the reflection of a secondary ray, only the reflected main ray is considered as the secondary ray. 

[Fortsetzung auf der ndchsten Seite] 



WO 02/33455 A2 IIUUIIIIIHIMHI 



— Erfindererklarung (Regel 4.17 Ziffer iv) nurjur US Zur Erkltirung der Zweibuchstaben-Codes und der anderen 

Abkurzungen wird auf die Erkldrungen ("Guidance Notes on 
VerdfTentlicbt: Codes and Abbreviations") am Anfang jeder reguldren Ausgabe 

— ohne international Recherchenbericht und erneut zu ^ er PCT-Gazette verwiesen. 
verqffentlichen nach Erhalt des Berichts 



(57) ZusammenfessuDg: Verfahren zur rechnerischen Bestimmung optischer Eigenschaften eines Kanalwellenleiters, bei dem fur 
einen in eine Eintrittsflache einfallenden Strahl die Intensitatsverteilung liber eine Austrittsflache durch aneinander anschlieBende, 
sich mdglicherweise aufteilende Teilstrahlen bestimmt wird. Aus einem auf die Mantelflache fallenden Teilstrahl wird durch ei- 
nen vorgegebenen Algorithmus ein reflektierter Hauptstrahl und, sofem nach dem folgenden Ablauf erforderlich, eine Anzahl von 
Streustrahlen hflherer Ordnung bestimmt, die ihrerseits TeilstrahJen darstellen und rekursrv weiterverfolgt werden. Jeder Teilstrahl 
ist entweder ein Primarstrahl oder ein Sekundarstrahl. Wenn der einfallende Teilstrahl ein Primarstrahl ist, ist der reflektierte Haupt- 
strahl wieder ein Primarstrahl und die Streustrahlen sind SekundMrstrahlen. Bei der Reflexion eines Sekundarstrahls wird nur der 
reflektierte Hauptstrahl als Sekundarstrahl berucksichtigt. 
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Beschreibung 

Strahlverfolgung in hoch multimodalen Kanalwellenleitern 



Tftchni arches Gehiet 

5 Die Erfindung betrifft die. Strahlverfolgung in hoch multimo- 
dalen Kanalwellenleitern. 

Zur Berechnung der Strahlausbreitung in dielektrischen Kanal- 
wellenleitern, insbesondere Lichtleitern, stehen bislang wel- 

10 lenoptische Analysemethoden wie die Methode der finiten Ele- 
mente (FEM) oder die 'Beam Propagation Method 1 (BPM) zur Ver- 
fttgung. Diese sind jedoch nur dann effizient einsetzbar, wenn 
nur eine oder wenige Moden zu berUcksichtigen sind und der 
Querschnitt der Wellenleiter, bezogen auf die optische Wel- 

15 lenlange, nicht allzu grofl ist. 

Fttr multimodale Stufenindex-Wellenleiter, bei denen der Quer- 
schnitt wesentlich grdfier als die Wellenl&nge der verwendeten 
Strahlung ist, sind hingegen Strahlverfolgungen auf der Basis 
20 geometrischer Optik effizient mQglich. 

Hierbei wird (in der Simulation) ein einzelner Strahl vorge- 
gebener Richtung und Polarisation in den Wellenleiter einge- 
koppelt. Dieser tritt entweder direkt am Ende des Wellenlei- 
25 ters aus oder wird an der Wand des optischen Kanals, d.h. der 
Grenzfl&che des Indexsprungs, gebrochen. 

Dabei teilt sich die eintreffende Strahlung in Komponenten 
auf: einen ref lektierten Hauptstrahl, einen transmittierten 
30 Hauptstrahl, eine Anzahl von ref lektierten Streustrahlen und 
eine Anzahl von transmittierten Streustrahlen. Die transmit- 
tierten Anteile sind ftir die weitere Strahlverfolgung ohne 
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Belang; lediglich ihr Energieanteil geht den ref lektierten 
Strahlen verloren. 

Eine einfache Simulation berticksichtigt lediglich den reflek- 
5 tierten Hauptstrahl, d.h. die Reflexion nullter Ordnung, und 
verfolgt dessen weitere Reflexion bis zur Austrittsf lache . 
Damit kann durch individuelle Verfolgung einer grOfieren An- 
zahl von Strahlen, die der Charakteristik des Senders ent- 
sprechen, das StrahlenbUndel am Ausgang des Wellenleiters be- 
10 stimmt werden. Diese Methode ist ausreichend/ wenn die Refle- 
xionen nahezu ideal sind, weil die Wandung sehr glatt ist, 

Ist jedoch die Wandung nicht glatt, so fiihrt die genannte Me- 
thode zu Ergebnissen, die entsprechenden Messungen nur 
15 schlecht entsprechen. Die Berttcksichtigung der ref lektierten 
Streustrahlung jedoch fiihrt zu einem Rechenaufwand, der expo- 
nentiell mit der Anzahl der Reflexionen steigt. 

Aufgabe der Erfindung ist es daher, eine Methode anzugeben, 
20 bei der auch die Streustrahlung berticksichtigt wird, dessen 
Rechenaufwand aber allenfalls quadratisch mit der Anzahl der 
Reflexionen steigt. 

Die Erfindung lost diese Aufgabe, indem sie zwischen Primar- 
25 und Sekundarstrahlen unterscheidet und bei Sekundarstrahlen 
nur noch Reflexionen nullter Ordnung berticksichtigt . 

Es handelt sich urn ein Verfahren zur rechnerischen Bestimmung 
optischer Eigenschaf ten eines Kanalwellenleiters, bei dem flir 

30 einen in eine Eintrittsf lache einfallenden Strahl die Inten- 
sitatsverteilung Uber eine Austrittsf lache durch aneinander 
anschlieftende, sich moglicherweise aufteilende Teilstrahlen 
bestimmt wird. Aus einem auf die MantelflSche fallenden Teil- 
strahl wird durch einen vorgegebenen Algorithmus eine reflek- 

35 tierter Hauptstrahl und, sofern nach dem folgenden Ablauf er- 
forderlich, eine Anzahl von Streustrahlen h5herer Ordnung be- 
stimmt, die ihrerseits Teilstrahlen darstellen und rekursiv 
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weiterverfolgt werden. Jeder Teilstrahl ist entweder ein Pri- 
marstrahl oder ein Sekundarstrahl . Wenn der einfallende Teil- 
strahl ein PrimSrstrahl ist, ist der reflektierte Hauptstrahl 
wieder ein Primarstrahl und die Streustrahlen sind Sekundar- 
5 strahlen. Bei der Reflexion eines Sekundarstrahls wird nur 
der reflektierte Hauptstrahl als Sekundarstrahl berticksich- 
tigt. 

TCiirghgsrhrgihimg der Efiichnunaen 
Es zeigen 

10 Fig. 1 eine zweidimensionale Darstellung eines Modells 
der Strahlverfolgung bei einem Kanalwellenleiter 
mit rauher Oberflache, 

Fig. 2 eine zweidimensionale Darstellung, reduziert auf 
die in der Erfindung verwendeten Teilstrahlen. 



15 BPRffhrPihing piner Ansf ilhrnngs f orm der FIrfinduna 

Fig. 1 zeigt eine zweidimensionale Darstellung zur Erlaute- 
rung des bisher verwendeten Verfahrens, wie es in dem Artikel 
"An Efficient Monte Carlo Based Ray Tracing Technique for the 
Characterization of Highly Multimode Dielectric Waveguides 
20 with Rough Surfaces" von Th. Bierhoff et al., Proc. 30th Eu- 
ropean Microwave Conference, Vol.1 pp 378-381, Paris, 3rd-5th 
October 2000, beschrieben ist. 

In Fig. 1 ist ein Langsschnitt durch einen Kanalwellenleiter 
25 20 skizziert, der einen Kern 21 und eine UmhUllung 22, auch 
als 'cladding 1 bezeichnet, hat. Der (nicht gezeigte) Que'r- 
schnitt kann rund oder rechteckig sein. Andere Querschnitte 
sind gleichfalls verwendbar, aber weniger praktisch in der 
Herstellung. Es handelt sich urn einen Stufenindex- 
30 Wellenleiter, weil ein relativ scharfer Obergang zwischen dem 
Kern 21 und der UmhUllung 22 vprhanden ist, so dafi eine Stufe 
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im Br echungs index entsteht, wobei bekanntermaflen der Bre- 
chungsindex des Kerns 21 grofier ist als der der Umhtillung 22. 
Dadurch entstehen Grenzfiachen 23a und 23b, an denen ein 
Strahl im Kern 21, der auf die Grenzflache trifft, reflek- 
5 tiert wird. Die Grenzfiachen sind hier gewellt dargestellt, 
urn anzudeuten, dail eine Oberfiachenrauhigkeit von Bedeutung 
ist; in dem behandelten Verfahren werden sie jedoch zunachst 
als Ebenen oder sonstige einfach, insbesondere analytisch zu 
beschreibende Fiachen angenommen, 

10 

Der Kanalwellenleiter hat eine Eintrittsf lache 24 und eine 
Austrittsfiache 25. In die Eintrittsf lache 24 failt ein 
(Licht-) Strahl 26, der an der Eintrittsf lache 24 gebrochen 
wird und als Eingangsstrahl 27 fortgesetzt wird. Ob und wie 
15 grofi die Brechung ist, hangt von der Art und dem Medium der 
Ankopplung ab; ftir das weitere Vorgehen ist lediglich der 
Eingangsstrahl 27 von Bedeutung. Dieser erreicht die Grenz- 
flache 23b an einem ersten Ref lexionspunkt 30 und wird dort 
gebrochen und reflektiert. 

20 

Die Berechnung eines solchen Strahlengangs mit den Mitteln 
der geometrischen Optik ist weithin bekannt. Hat der Kanal- 
wellenleiter rechteckigen Querschnitt, dann sind die Grenz- 
fiachen 23a, 23b wiederum rechteckig, so dafi bei gegebener 

25 Gerade fUr den Eingangsstrahl . 27 der Schnittpunkt der Geraden 
mit einer der Grenzfiachen, hier der Grenzflache 23b, mit den 
bekannten Methoden der analytischen Geometrie bestimmt werden 
kann. Ist der Querschnitt rund, so ist die Grenzflache 23a, 
23b ein Zylinder, zu dem der Schnittpunkt auch noch einfach 

30 berechnet werden kann. Gegebenenf alls mufi bei ungewohnlichen 
oder variablen Querschnitten auf iterative numerische Metho- 
den zuruckgegrif fen werden. 

Die Fortsetzung des am ersten Ref lexionspunkt 30 eintreffen- 
35 den Eingangsstrahls 30 richtet sich nach den Parametern der 
Simulation. Danach wird die Intensitat und die Richtung des 
direkt ref lektierten Strahls 31a bestimmt. Die Richtung rich- 
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tet sich nach der Tangente der idealisierten Grenzfiache 23b, 
die als rauh angesehen wird. Ferner werden diffuse Streu- 
strahlen bestimmt, die in Fig. 2 durch einen einzigen Streu- 
strahl 31b reprasentiert werden. Hierbei wird ein Raumwinkel, 
5 in dem Streustrahlung kontinuierlich vorhanden ist, zu einem 
einzigen Strahl in der Modellierung zusaramengef afit . Der 
Streustrahl 31b reprasentiert damit die Streustrahlung in ei- 
nem gegebenen Raumwinkel, die durch eine ebene Welle mit de- 
finierter Polarisation transportiert wird; fttr benachbarte 

10 Raumwinkel werden entsprechend weitere Streustrahlen be- 
stimmt, die auf gleiche Art weiter behandelt werden, aber der 
Ubersichtlichkeit halber nicht dargestellt sind. Gleichfalls 
angedeutet durch Pfeile sind ein gebrochener Primarstrahl und 
ein gebrochener Streustrahl. Deren weiterer Weg ist ohne Be- 

15 deutung; lediglich ist f estzuhalten, dafi die Summe aller In- 
tensitaten gleich der des eintref fenden Eingangsstrahls 27 
sein mufi. 

Der direkt reflektierte Hauptstrahl 31a trifft an der Stelle 
20 32 wieder auf eine Grenzfiache, hier die Grenzfiache 23a, und 
wird dort zu einem ref lektierten Hauptstrahl 33a und mehreren 
reflektierten Nebenstrahlen 33b, von denen wieder nur einer 
dargestellt ist. Die transmittierten Strahlen werden wieder 
nur ihrem Anteil nach berticksichtigt . 

25 

Entsprechend trifft der reflektierte Streustrahl an der Stel- 
le 36 auf die Grenzfiache 23b und ergibt einen Hauptstrahl 
und mehrere Streustrahlen. 

30 Ersichtlich ist dieses Verhalten von der Struktur durch re- 
kursive Verfahren relativ einfach programmierbar . Es ist je- 
doch auch sofort ersichtlich, daB insbesondere bei mehreren 
Streustrahlen die Anzahl der Strahlen exponentiell steigt. 
Dabei werden Strahlen, deren Intensitat unter einer vorgege- 

35 benen Schwelle liegen, natUrlich nicht weiter berticksichtigt. 
Man muB jedoch feststellen, daB entweder die Schwelle so hoch 
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liegt, dafi die Erg.ebnisse nicht gut genug sind, Oder zuviel 
Rechenzeit benStigt wird. 

In dem genannten Artikel wird daher vorgeschlagen, die Kom- 
5 plexitat durch Anwendung der Monte-Carlo-Methode zu reduzie- 
ren. Diese allgemein bekannte Methode wahlt durch einen Zu- 
fallsmechanismus einen der Streustrahlen aus und ignoriert 
die anderen. Es ergab sich zwar eine Reduktion der Rechen- 
zeit, aber insgesamt noch keine einfache, schnelle und prob- 
10 .lemlos benutzbare Lbsung. 

Durch die Erfindung kann auf die Monte-Carlo-Methode verzich- 
tet werden und ein wesentlich besseres Ergebnis in Relation 
zum Rechenaufwand erreicht werden. Dabei werden alle nach dem 

15 Modell der Reflexion entstehenden Streustrahlen weiterver- 
folgt und nicht nur einige zufailig ausgewahlte. FUr diese 
Streustrahlen jedoch wird bei nachf olgenden Reflexionen nur 
noch der Hauptstrahl berticksichtigt . Dabei wird jedesmal das 
Modell der Reflexion berechnet, urn die Intensitat des Haupt- 

20 strahls zu ermitteln; es wird also nicht unbedingt ein einfa- 
ches Spiegelmodell fiir diese Nebenstrahlen eingesetzt, obwohl 
dies auch mttglich ist. 

Dies wird erreicht, indem zu jedem zu verfolgenden Strahl zu- 
25 satzlich ein Generationsindex.mitgefilhrt wird, der angibt, ob 
es sich urn einen Primarstrahl oder einen SekundcLrstrahl han- 
delt. Bevorzugt wird hierzu eine ganze Zahl benutzt, die Null 
ftir einen Primarstrahl und grolier Null ftir einen Sekundar- 
strahl ist. 

30 

Dies ist in Fig. 2 skizziert. Die transmittierten Anteile 
wurden ohnehin nicht mehr dargestellt. Aber auch die von dem 
sekundaren Streustrahl 31b bewirkten sekundaren Streustrahlen 
werden gar nicht mehr weiterverfolgt . Hingegen wird der se- 
35 kundare Teilstrahl 33b als von einem primaren Teilstrahl 31a 
bei der Reflexion 32 ausgehend weiterhin bertlcksichtigt. 
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Bei der einfachen, bislang beschriebenen AusfUhrung wird le- 
diglich ein binarer Generationsindex benutzt bzw. nur zwi- 
schen Printer- und Sekundarstrahlen unterschieden. Es ist je- 
doch auch moglich, die Simulation zu verbessern, indem der . 
5 Generationsindex als ganze Zahl ausgenutzt wird. Dabei bleibt 
der Generationsindex bei dem Hauptstrahl einer Reflexion, 
d.h. der Reflexion nullter Ordnung, gleich. Bei Streustrah- 
len, d.h, Reflexionen httherer Ordnung, wird die Ordnungsnum- 
mer zu dem Generationsindex addiert. Ob ein Strahl weiterver- 
10 folgt wird, wird dann Uber einen Schwellwert fUr den Genera- 
tionsindex entschieden. In dem oben genannten einfachen Bei- 
spiel werden Uberhaupt nur Reflexionen erster Ordnung berUck- 
sichtigt, und der Schwellwert ist Eins. 

15 Durch den Gewinn an Effizienz konnte bei der Reflexion auch 
die Polarisation berUcksichtigt werden, was bei der aufge- 
ftlhrten Monte-Carlo-Methode ignoriert werden mulite. 
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Patentanspriiche 

1. Verfahren zur rechnerischen Bestimmung optischer Eigen- 
schaften eines Kanalwellenleiters (20) mit einer Mantel- 
fiache, wobei ftir einen in eine Eintrittsf l^che (24) ein- 

5 fallenden Strahl (27) die resultierende Strahlverteilung 
iiber eine Austrittsf lache (25) bestiimnt wird, mit den 
Merkmalen: 

- zu einem auf die Mantelflache (23a, 23b) fallenden 
Teilstrahl (31a, 31b, 33a, 33b, 35a, 35b) wird ein re- 

10 flektierter Hauptstrahl (31a, 33a, 35a) der Ordnungs- 

zahl Null und, sofern nach dem folgenden Ablauf erfor- 
derlich, eine Anzahl von Streustrahlen (31b, 33b, 35b) 
mit einer Ordnungszahl grSSer Null bestiramt, die ih- 
rerseits Teilstrahlen darstellen und rekursiv weiter- 

15 verfolgt werden, wobei der einfallende Strahl (27) der 

erste Teilstrahl ist, 
dadurch gekennzeichnet , daS 

- jedem Teilstrahl ein Generations index zugeordnet ist, 

- der einfallende Strahl den Generationsindex Null hat, 

20 - bei der Reflexion der Generationsindex eines reflek- 

tierten Teilstrahls gleich dem um die Ordnungszahl der 
Reflexion erhohten Generationsindex des einfallenden 
Teilstrahls ist, 

- nur solche Teilstrahlen weiterverf olgt werden, deren 
25 Generationsindex eine vorgegebene Schranke nicht iiber- 

schreitet . 

2. Verfahren nach Anspruch 1, wobei die vorgegebene Schranke 
gleich Eins'ist und bei der Berechnung der Reflexion al- 
lenfalls nur Streustrahlen der Ordnungszahl Eins anfal- 

30 len. 
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